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Abstract 

We present the first phenomenological study of the masses of orbitally excited 
baryons in large N c QCD. Restricting here to the nonstrange sector of the 
t = 1 baryons, the 1/N C expansion is used to order and select a basis of 
effective operators that spans the nine observables (seven masses and two 
mixing angles). Fits are performed using subsets of the complete set of nine 
operators, including corrections up to 0(1/N C ) where leading order is JV~. 
This study shows that the 1/N C expansion provides an excellent framework 
for analyzing the mass spectrum, and uncovers a new hierarchy of operator 
contributions. 



I. INTRODUCTION 



It appears that QCD admits a useful and elegant expansion in powers of 1/N C , where 
N c is the number of colors 0. There are explicit rules that determine the order in N c of 
any given Feynman diagram or matrix elements of any given operator. One can thus isolate 
first the leading contributions to the observable under consideration and then systematically 
include contributions that are proportional to higher powers of 1/N C . 

One may question whether 1/N C = 1/3 is small enough to provide a valid phenomeno- 
logical expansion parameter. Experience suggests that it is. For the ground state baryons, 
the large N c approach has been used successfully to study SU(6) spin-flavor symmetry |2]-[|. 
masses , magnetic moments H,|8|JToHl2| , and axial current matrix elements [pj^|,|8|,[l2[ . 



The next natural step in this progression is studies of excited baryons, in particular of 
the mixed symmetry, negative parity 70-plet of SU(6). There have been studies of the strong 
13| , |T4|| and radiative [15[ decays of these states and of the structure of axial operator matrix 



elements [ 16| . Progress [|17| along the lines of the present work has already been made with 
the first non-leading mass operators. 

In this paper, we study the relative order 1/N C and 1/N% corrections to the masses of 
the nonstrange members of the 70-plet. Qualitatively, we find that the large N c limit is 
accurate for the 70-plet, in the sense that all the mass operators — if one wishes, all the mass 
terms in an effective Hamiltonian — give contributions at or below the level estimated from 
large N c considerations. None are larger. Quantitative detail is added to this statement in 
Section III. 

Before proceeding, we make some comments on the nature of the 70-plet and the nomen- 
clature we use. We describe the state as a symmetrized "core" of (N c — 1) quarks in the 
ground state plus one excited quark in a relative P state. The wave function is antisym- 
metric in color and symmetric in SU(6) x 0(3), where SU(6) is the spin-flavor symmetry 
and 0(3) is the rotation group. We use SU(6) to classify the large N c baryon states and 
the transformation properties of the operators; however, we need not assume that SU(6) is 
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an exact symmetry. In fact, while the leading order contribution to the masses is of 0(N C ), 
SU(6) is broken at 0(N®) for our states. 

One can analyze the masses in the 70-plet by expressing the effective Hamiltonian as 
a sum of operators, one of which is proportional to the identity and the rest are products 
of SU(6) x 0(3) generators times numerical coefficients. Each operator contributes to the 
mass at a definite order in N c , which is determined by rules delineated in Sec. ||. 

The operator analysis for the excited baryons involves more distinguishable generators 
than the ground state baryons, because of the orbital angular momentum and because one 
quark is singled out. This leads to a much larger collection of mass operators. A list is given 
in the next section, which shows that there is one operator of 0(N C ), two operators of 0(1), 
and many of 0(1/N C ) or smaller. 

Regarding the states, nonstrange mixed symmetry SU(6) states with one quark singled 
out have total quark spin and isospin related by S = I or / ± 1, with each of S and / in the 
range 1/2 to N c /2. (There is one exception: There are no doubly maximal mixed symmetry 
S = I = N c /2 states.) Thinking of the states as a core of N c — 1 ground state quarks with 
spin S c and isospin I c combined with an excited quark with angular momentum i = 1 leads 
to writing the SU(6) x 0(3) states, with the help of Clebsch-Gordan coefficients, as 
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Here, the m's are angular momentum projections, the a's are isospin projections. Note that 
S c = I c since the (nonstrange) core is symmetric in SU(6) indices, and we have written 
I c — I + 77/2, where rj = ±1. Note also that p = S — I = ±1,0. States with strangeness 
are defined analogously, except that SU(3) Clebsch-Gordan coefficients appear in that case. 
For reasons of simplicity, we have restricted to nonstrange baryons in this work. 
For S = I ± 1, the notation simplifies since c± t ± = 1 and c± iT = 0. For S = I, 
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c 0+ = + , 



S(N C + 2(S + 1)) ^ ^ 



(S + 1)(N C -2S)_ 



\ N C (2S + 1) u "~ \ N C (2S + 1) ■ 

(The orthogonal combination gives the totally symmetric SU(6) state.) The explicit form 
of the state allows us to calculate analytically the matrix elements for any operator for 
arbitrary iV c . 

For the physical case of N c = 3, the above expressions admit the 7 nonstrange states of 
the 70-plet. Strictly speaking, the label 70-plet refers only to the mixed symmetry baryons 
appearing at N c = 3; for large N c , the representations tend to be much larger. In the 70-plet, 
the nonstrange states consist of two isospin-3/2 states, A1/2 and A3/2, and five isospin-1/2 
states, N1/2, N[/ 2 , N 3 / 2 , NL 2 , and NL 2 . The subscript indicates total baryon spin; unprimed 
states states have quark spin 1/2 and primed states have quark spin 3/2. In partial wave 
notation, the 2 deltas and 5 nucleons are labeled S31, -D33, Su (twice), (twice), and Di 5 , 
respectively. 

Section || explains the operator analysis and gives matrix elements of a complete set of 
operators for the baryon states. Section |l| contains our analysis of physical masses and 



mixing angles. Closing comments are made in Sec. IV. 



II. OPERATORS AND MATRIX ELEMENTS 

Since the physically observed baryons are assigned to irreducible representations of the 
symmetry group SU(6) x 0(3), it is natural to write all possible mass operators in terms 
of the generators of this group. 0(3) is the group of spatial rotations and is generated by 
the orbital angular momentum operator £\ while the spin-flavor group SU(6) has spin S\ 
flavor T a , and combined spin-flavor G ia generators. In the two-flavor case, these operators 
are defined by 
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G m = g f I y ® T J g, (2.1) 

where <r* and r a are the usual Pauli matrices. The field operators q, q\ which we term 
"quarks," are not the dynamical quarks, but rather eigenstates of the spin-flavor group such 
that an appropriately symmetrized collection of N c of them have the quantum numbers of 
the physical baryons. The collection of operators constructed with these fields completely 
spans all possible physical mass operators. Only if the quarks are heavy can the fields q be 
identified with the dynamical valence quarks. 

The goal of the large N c analysis is to organize operators by their effects on a given 
observable (in this case, the masses) in a systematic expansion in powers of N c . Factors 
of N c originate either as coefficients of operators in the Hamiltonian, or through matrix 
elements of those operators. For example, the unit operator IL contributes at O(N^), since 
each quark contributes coherently in the matrix element. The spin of the baryon S 2 is known 
to contribute to the masses at 0(1/N C ) 0, because the matrix elements of 5** are of 0(N®) 
for baryons that have spins of order unity as N c — > oo. Similarly, matrix elements of T a are 
0(N®) in the two-flavor case since the baryons considered have isospin of 0(N®), but the 
operator G m has matrix elements on this subset of states of O(N^). This means that the 
contributions of the N c quarks add incoherently in matrix elements of the operator S' 1 or T a 
but coherently for G m . Note that matrix elements of a given operator are not necessarily 
homogeneous in N c ; for example, values such as N c + 3 can occur. 

In this work, the generators SI, T c a , G* a are reserved to mean those acting upon the 
core, while separate SU(6) generators s % , t a , and g m are defined for the single excited-quark 
system. Including the operator t completes the list of building blocks for the necessary mass 
operators. Strictly speaking, the naive symmetry group for this operator basis is SU(6) x 
SU(6) x 0(3), although actually only the diagonal subgroup SU(6) x 0(3) truly acts on 
the baryon states. Given this enlarged notation, it is possible to construct a large number 
of operators; however, many are linearly dependent and may be discarded. 

Since the excited system consists of only one quark, at most one generator from among 
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s, t, or g appears in any operator. A similar but stronger statement may be made for the 
core: Since we ultimately perform phenomenological analysis on cores with two quarks, at 
most two of the set {S c , T c , G c } are needed. However, operator reduction rules exist || that 
significantly reduce the number of core operators that must be considered. Finally, since 
we are interested here in I = 1 baryons, only operator combinations of A£ = 0,1,2 need be 
considered. Hence, only up to two factors of t are required. Indeed, when two ts appear, 
it is convenient to use l\j , the rank two tensor combination of t and i 1 given by 

4 2) = ^{4^-} - f (2.2) 

The explicit power of N c for a given operator is determined by using the usual large 
N c counting of 1/ \fW c for each quark-quark-gluon coupling, and considering the minimal 
number of exchanged gluons necessary to generate a given operator. To be specific, we 
decompose an n-body operator O as follows: 

71-1 

O = X* J] *i, (2-3) 

i=l 

where X* represents all operators acting on the excited quark, including factors of t, and 
each Xi represents an SU(6) generator acting on the core. The physical realization of such 
an operator requires exchanging a minimum of n — 1 gluons between different quarks, leading 
to a suppression of 1/iV" -1 , which we henceforth include in the definition of the operators. 
If X* = 1, the result is maintained as written if n rather than n — 1 operators act on the 
core. 

One then considers each core operator Xj to determine whether its matrix elements 
are coherent for the baryon states under consideration; a factor of N c is included for each 
coherent operator. Thus, the full large N c counting of the matrix element is 0(X c 1_ra+m ), 
where m is the number of coherent Xj. In the two-flavor case, only G* a has coherent matrix 
elements. The order of the matrix element thus obtained determines whether or not we 
retain the operator in computing a given process to a desired order in the 1/N C expansion. 

With this counting, one finds 22 potentially independent, time- reversal even, isosinglet 
operators: One (1) with matrix element of O(N^), two at 0(N®), ten at O(l/N ), and nine 
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at 0(1/N%) or higher. This counting does not fully take into account numerous relations 
between the matrix elements of the operators evaluated on the nonstrange I = 1 baryons. 
One reduction that has been included uses the observation that, for the nonstrange mixed 
symmetry states, 

N c W cl 16iV c bJ 2N* V ; 

This operator naively produces 0(N®) matrix elements, but both the 0(N®) and 0(1/N C ) 
parts are the same for all baryons in the multiplet, and thus may be absorbed into matrix 
elements of 11, with the remainder being demoted to 0(1/N%). Similarly, both (is) and 
(£tG c ) are O(iV c ), but it may be observed that {Is + UtG c /(N c + 1)) is 0(1/N C ), so only 
(is) truly represents an independent 0(N®) operator. The full set of operator reductions 
will be presented in a future publication . 

In any case, there are only 9 observables (masses and mixing angles) in the nonstrange 
1 = 1 system, and so only 9 independent operators are required. An independent basis 
is presented in Table Itt All 3 occurring up to 0(N®) (first considered in Ref. [0), and a 
selection of 6 at 0(1/N C ) whose matrix elements, when combined with the first 3, are seen to 
be independent for N c = 3. With these operators denoted by 0%, O2, ■ ■ ■ , Og, respectively, 
the nine independent mass matrix elements are given by 

X'v(a), (j = 1...9). (2.5) 
i=i 

The coefficients q are independent of N c at leading order, given our choice of operator nor- 
malization. These operator coefficients encapsulate all unknown strong interaction physics 
unspecified by the large N c spin-flavor analysis. In Table | we present the matrix elements 
(Oi)j. Explicit spin and flavor indices are suppressed when their contraction is unambiguous. 



III. RESULTS 



In addition to the nonstrange mixed symmetry states defined in Sec. |, two mixing angles 
are necessary to specify the 5=1/2 and 5 = 3/2 nucleon mass eigenstates. We define 
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as in Ref. [[Uj . The mass eigenvalues and mixing angles can be expressed in terms of the 
coefficients of the operators presented in the previous section. 

Since we have found an operator basis that completely spans the 9-dimensional space 
of observables, we can solve for the Cj given the experimental data. For each baryon mass, 
we assume that the central value corresponds to the midpoint of the mass range quoted in 
the Review of Particle Properties ||19|| ; we take the one standard deviation error as half of 
the stated range. To determine the off-diagonal mass matrix elements, we use the mixing 
angles extracted from the analysis of strong decays given in Ref. 9^i = 0.61 ± 0.09 and 
9n3 = 3.04 ± 0.15. These values are consistent with those obtained in [HjJ] from radiative 



decays. Solving for the operator coefficients, we obtain the values shown in Table [TT|. 

Naively, one expects the Q to be of comparable size. Using the value of c\ as a point of 
comparison, it is clear that there are no operators with anomalously large coefficients. Thus, 
we find no conflict with the naive 1/N C power counting rules. It is interesting that a number 
of the operators appear to be unimportant in describing the experimental data (presumably 
due to the underlying dynamics). For example, of the two operators that contribute to 
the masses at 0(1), the operator O2 = is has a coefficient which is suppressed relative to 
(9 3 = &( 2 >gG c /N c by more than factor of 10; in effect, this operator is no more important 
than a typical 0(1/ N 2 ) correction. Of the operators that contribute to the masses at order 
1/N C , only the operator Oq = S^/N c contributes as much as one would expect, with the next 
largest corrections coming from the operators C 4 = is + 4itG c /(N c + 1) and C 5 = iS c /N c . 

Using these observations, we can attempt to fit the data using judiciously chosen subsets 
of the original 9 operators. We fit to the seven mass eigenvalues as well as the two mixing 



angles 6ni and 9 m- The operator set we consider first are Ox, 0%, and O3, these yield mass 
predictions accurate to order 1 in the 1/N C expansion, and thus present the first nontrivial 
spin-flavor symmetry-breaking corrections. We show the result of this fit in Table |TJ. In 
short, the lowest order operators fail in reproducing the experimental data. Notably, the 
mixing angles are far off the mark, and the J = 1/2 A state is predicted to be heavier 
than the J = 3/2 A state. The difficulty in obtaining a good fit from a leading order 
analysis is an outcome that perhaps could have been anticipated: Naively, one might expect 
that the lowest nontrivial mass corrections are roughly a factor of N c = 3 smaller than the 
mean baryon mass, or approximately 500 MeV. The largest splitting within our set of seven 
baryons is ~ 180 MeV, for example, in the case of the N(1700)-N(1520) mass difference. 
Thus, we might have concluded a priori that 1/N C corrections are necessary in order to 
reproduce the detailed features of the mass spectrum. 

In the remaining fits, we include 1/N C corrections. The 6 parameter fit shown in Table p"V| 
includes all the subleading operators that appear to be significant in Table [TI]; the operators 
included are 0\ through Oq. The resulting fit is in extremely good agreement with the 
experimental data, with no predicted mass more than 0.4 standard deviations from the 
corresponding experimental central value, and a x 2 P er degree of freedom of 0.1. 

More strikingly, Table [ll| implies that we can discard additional operators and still obtain 
a reasonable fit. Notice that the smallness of the coefficients c-i, C4, and C5 renders the 
corresponding operators numerically unimportant, and thus they can be neglected if we are 
only interested in working to order 1/N C . In Table [V| we give a fit retaining the remaining 
three operators, Ox, O3, and Oq. The \ 2 P er degree of freedom for this fit is 1.87, which is not 
bad considering that have only included two nontrivial operators. Notice that the particular 
operator choice for this fit leads to a degeneracy between the A1/2 and A 3 / 2 which is lifted 
by the corrections that we have discarded. We do not display the fits corresponding to all 
possible choices for the subleading operators. It suffices to point out that these additional 
fits interpolate between the 6 and 3 parameter fits that we have presented in Tables [TV] 
and [V|. One may also fit to the mass eigenvalues and predict the mixing angles. These fits 
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are not qualitatively different from the ones given here, and will be presented in a longer 
publication |18[. The crucial observation is that the sub leading operator C 6 = S 2 /N c is the 
most significant ingredient in taking us from the poor fit shown in Table [TIT] to the good fits 
in Tables [V] and [V]. 



IV. CONCLUSIONS 

The value of the large N c approach to baryon phenomenology is that it provides an 
organizing principle in constructing the baryon effective field theory. Studies of the excited 



baryon mass spectrum in the formative days of SU(6) found numerous operators |20[], but 
in that period there was no organizing principle available to select among them. Beginning 
with a complete operator basis that spans the space of any desired set of observables, large 
N c power counting rules tell us which operators may be discarded if we wish to obtain 
predictions to a desired level of accuracy. In this sense, our results for the masses and 
mixing angles of the nonstrange I = 1 baryons are completely consistent with the large N c 
picture. We find no operator with a coefficient that is larger than what one would expect 
from the naive large N c power counting rules. 

More interesting, however, is that at any given order in 1/N C , only some of the operators 
are of phenomenological relevance. The 0(N®) coefficients Cj that we have defined in our 
effective theory parametrize the long-distance physics that we cannot calculate. If we had 
found that these coefficients were of comparable size, we might have concluded that large 
N c counting arguments alone are sufficient to explain the detailed features of the mass 
spectrum. Quite the contrary, we find that only a few of our original set of operators are 
needed to reproduce the experimental data, most notably, the operators S 2 and £^ 2 'gG c . 
It is tempting to speculate that the importance of the operator S 2 can be understood by 
considering the explicit nonrelativistic reduction of a one-gluon exchange interaction; the 
second operator, however, has a nontrivial flavor structure that does not correspond to the 
usual tensor interaction £( 2 >sS c that one derives in this approach [ETJ . A fit analogous to 
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that of Table [V] replacing £( 2 'gG c /N c with ^ 2S} sS c /N c increases the x 2 P er degree of freedom 
from 1.87 to 2.46. Why the underlying dynamics should prefer these operators is a much 
more difficult question which goes beyond what can be addressed in the effective field theory 
approach. 
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TABLE I. Matrix elements (Oi)j of 9 operators, labeled as 0±, O2, ■ ■ ■ , Og, respectively, that 
are linearly independent for N c = 3. The third and sixth rows correspond to off-diagonal matrix 
elements. 
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TABLE II. Operator coefficients in GeV, assuming the complete set of Table Q. The vertical 
divisions separate operators whose contributions to the baryon masses are of orders ./V,, 1 , N® and 
N~ , respectively. 



Parameters (GeV): c% = 0.542 ± 0.002, c 2 = 0.093 ± 0.008, c 3 = -0.335 ± 0.039 
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TABLE III. Three parameter fit using operators Ci,2,3j giving x 2 /d.o.f. = 23.33/6 = 3.89. The 



operators included formally yield the lowest order nontrivial contributions to the masses in the 
1/N C expansion. 
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Parameters (GeV): a = 0.468 ± 0.005, c 2 = -0.032 ± 0.045, c 3 = 0.327 ± 0.093 



c 4 = 0.081 ± 0.027, c 5 = 0.071 ± 0.042, c 6 = 0.413 ± 0.044 
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TABLE IV. Six parameter fit using operators Oi t .. 


,6, giving x 2 /d.o.f. 


= 0.31/3 = 0.10. 


Parameters (GeV): a = 0.461 ± 0.005, c 3 = 0.360 ± 0.059, c 6 = 0.453 ± 0.030 




Fit (MeV) 


Exp. (MeV) 




Fit (MeV) 


Exp. (MeV) 


A(1700) 


1683 


1720 ± 50 


JV(1520) 


1530 


1523 ± 8 


A(1620) 


1683 


1645 ± 30 


JV(1535) 


1503 


1538 ± 18 


7V(1675) 


1673 


1678 ± 8 


&N1 


0.45 


0.61 ±0.09 


iV(1700) 


1725 


1700 ± 50 




3.04 


3.04 ±0.15 


7V(1650) 


1663 


1660 ± 20 









TABLE V. Three parameter fit using operators Oi, O3, and Oq, giving x 2 /A.o.i. = 



11.19/6 = 1.87. 
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